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A technique  la  present'd  .'or  obtaining  confidence  Intervals  and 
confidence  band*  for  a one-dinenslonai  response ’ surface  which  Is  assumed 
to  be  a general-order  polynomial.  Specifically,  suppose  that  Its  an 
Input  parameter  for  an  experiment  and  that  estimates  for  an  output  g( A ) 
are  desired,  where  g la  an  (n-I) th-order  polynomial  over  an  Interval 
*„+!•  If  obaervatlons  are  made  at  n •Ji«*-fr.wC  values  X-  at, 

chen  under  certain  conditions  confidence  Intervals  may  be  obtained 
for  g(X),  a^EX  f.a^j,  and  a confidence  band  may  be  obtained  for  the  entire  , 
function  over  this  Interval.  The  basic  requirement  for  these  results  Is 
that  central  limit  theorem*  exist  which  permit  confidence  interval  estima- 
tion at  the  n points  X.  a^,  two  Illustrations  are  given  which  satisfy 
this  requirement,  including  an  application  of  the  regenerative  approach  In 
steady-state  simulation  analysis. 
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ABSTRACT 


A technique  is  presented  for  obtaining  confidence  intervals  and 
confidence  bands  for  a one-dimensional  response  surface  which  is  assumed 
to  ,be  a general-order  polynomial.  Specifically,  suppose  that  X is  an 
input  parameter  for  an  experiment  and  that  estimates  for  an  output  g(X) 
are  desired,  where  g is  an  (n-1) th-order  polynomial  over  an  interval 
a^<X<  a . If  observations  are  made  at  n distinct  values  X-  a-, 
lg<i_<n,  then  under  certain  conditions  confidence  intervals  may  be  obtained 
for  g(X),  an£X  <an+^,  and  a confidence  band  may  be  obtained  for  the  entire 
function  over  this  interval.  The  basic  requirement  for  these  results  is 
that  central  limit  theorems  exist  which  permit  confidence  interval  estima- 
tion at  the  n points  X = a^.  Two  illustrations  are  given  which  satisfy 
this  requirement,  including  an  application  of  the  regenerative  approach  in 
steady-state  simulation  analysis. 
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TECHNIQUES  OF  RESPONSE  SURFACE  ESTIMATION, 

WITH  APPLICATIONS  IN  COMPUTER  SIMULATION* 

by 

Michael  A.  Crane 

i . Introduction 

In  this  paper  a method  is  presented  for  one-dimensional  response 

surface  estimation  when  it  is  assumed  that  the  response  surface  has  a 

general-order  polynomial  form.  In  particular,  it  is  shown  that  if  the 

response  is  of  order  n - 1 and  if  observations  are  taken  at  n distinct 

values  for  the  input  or  "treatment"  variable,  than  under  general  conditions 

it  is  possible  to  obtain  large-sample  confidence  intervals  for  the  response 

at  any  value  of  the  input  variable.  Additionally,  large-sample  confidence 

bands  for  the  entire  response  function  may  be  obtained. 

Specifically,  suppose  that  A is  an  input  parameter  for  a statistical 

experiment  and  that  estimates  for  an  unknown  output  function  g(A)  are 

desired,  where  g is  an  (n-lt^)-order  polynomial  over  an  interval 

a„  < A < a ...  Suppose  further  that  observations  are  made  for  n distinct 
0 — — n+1. 

values  of  A,  namely,  A = a , , 1 < i < n with  an  < a < ...<  a < a . 

i — — 0—1  n—  n+1 

In  Section  2,  conditions  are  given  which  allow  one  to  compute  confidence 

intervals  for  g(A)  at  any  point  an  < A < a , . In  Section  3,  confi- 

U — — n+1 

dence  bands  for  the  function  g(A)  over  a.  < A < a are  derived.  The 

0 — — n+i 

basic  requirement  for  these  results  is  that  central  limit  theorems  exist 
which  permit  confidence  interval  estimation  at  the  n points  A = a^. 

, * 

Thi3  research  was  sponsored  under  Office  of  Naval  Research  Contract 
N00014-72-C-0086  (NR-047-106) . 
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The  methods  are  illustrated  in  Section  4,  where  two  examples  are  provided 
which  satisfy  this  requirement.  The  first  example  is  the  typical  case 
where  the  experiment  at  A = consists  of  a sample  of  independent  and 
identically  distributed  variates  from  a population  with  mean  g(a^).  In 
the  second  example,  the  obseivations  at  A = are  based  on  random 
tours  in  an  application  of  the  regenerative  approach  in  simulation  analysis; 


see  ■ [4,  5,  6,  7,  10] . 

The  problem  as  stated  above  has  been  studied  previously  in  the  concext 
of  linear  regression  analysis.  The  interval  estimation  problem  is  discussed, 
for  example,  in  Graybill  [8,  pp.  121-122],  and  confidence  bands. for  a general 
linear  model  are  obtained  by  Bowden  [2].  However,  the  results  presented  here 
are  co  be  distinguished  from  rhe  linear  regression  models  in  that 

1)  the  confidence  interval  and  confidence  band  derivations  do  not 
require  distributional  assumptions, 

2)  the.  variances  at  A = a^  are  unknown  and  no  relationship  among 
the  variances  is  assumed, 

3)  the  response  surface  is  a general  polynomial,  and 

4)  it  Is  not  essential  that  all  observations  be  independent,  only 
that  n independent  experiments  be  conducted  and  that  a cential 
limit  theorem  be  applicable  for  the  observations  in  each  experi- 


2 . Confidence  Interval  Estimation 

Suppose  that  A is  an  Input  parameter  or  "treatment"  variable  for  a 
statistical  experiment  and  that  we  wish  to  estimate  some  output  parameter 
g(A)  which  is  an  unknown  (n-1)  t*1-order  polynomial  function  of  A over  an  in- 
terval 3q  £ A < an+-^«  i-e*. 


g(A)  - Zl  + x2  A ♦ ...  + zn  A , nQ  <A  < 


2 


where  ’-he  coefficients  z,  , z„,  z are  unknown.  Independent  experiments 

Z n 

reiulting  in  statistical  observations  are  made  at  n particular  parameter 
settings  X - a.,  1 < i £ n,  with  afl  < a^  < < ...  < a^  < a^. 

At  the  parameter  setting  A = a the  experiment  produces  observations 

resulting  in  sample  statistics  ^(a,,  N ) and  fi(a^,  N^) . for  example,  the 
N.  observations  might  he  independent  and  identically  distributed  sample  var- 
iates from  a population  with  mean,  g(a^),  and  ?(?^,  N^)  and  T)  (a^» 
might  be  the  sample  mean  and  sample  standard  deviation,  respectively.  Alterna 
tively,  the ' observat ions  migh : be  based  on  random  tours  in  an  application  of 
the  regenerative  approach  in  simulation  analysis.  See  Section  4 for  details. 
The  function  g may  be  expressed  in  terms  of  its  values  at  the  n 

points  a,,  a„ , ....  a , as  follows: 

1 Z n 


g(A)  = f (A ) 


g(ax) 


8(a„) 


• *0  - 1 - “n*l 


where 


V(A)  - (1  A ... 


x"-1) 


I “1 


n-l 

1 a x . . . 3^ 


1 a 


n-l 


n n 


Let  ^(A)  denote  the  component  of  y(A).  Note  that  the  matrix  to  be 

inverted  is  the  Vandermonde  matrix,  and  its  inverse  is  wel  1-known,  c f . [«,] 
pp.  124-125. 
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The  following  proposition  forms  a basis  for  a method  to  compute  confi- 
dence intervals  for  gO)  where  a < \ < a . . Let  a double  arrow  denote 

0 - - n+1 

convergence  in  distribution  and  let  N(p,<j)  denote  a normally  distributed 
random  variable  with  mean  p and  standard  deviation  <j  . 

PROPOSITION  1.  Suppose  there  esist  p(a^)  such  that  fj(a.,  N^)  > p(a  ) 


in  probability  as  N, 


d,  1 = 1,  2,  . ..,  n.  Suppose  further  that 


Ni[*(ai’  *V  ~ H(a.)  =*>  N(0,  1) 


•as  N.  — ^ 00 , i = 1,  2,  ....  n.  Then  for  a„  < A < a.,, 
— I U — — n't"  I 


[?(A,  n , ....  Nn)  - g ( A)  ] / v(A,  Nj Nn)=»N(0,  1) 


Ni— > CO  and  =*  [Ci  N^]  (Integer  part)  i * 2,  ...»  n,  where 


Cj > • • • » are  arbitrary  positive  constants. 


?( A , N. , ...»  N ) - V (A) 
in 


'f(alt  Mx) 


f(a  . N ) 


n’  Vj 


v(A,  # • • • * 


v - {5  'i(a)  «2<v 


Proof . From  the  conditions  of  the  proposition,  we  may  write 


[Ci  Nx  ]1/2  [9  (ai,  Nt)  - g^)]— * N(0,n(ai)),  i » 1,  ....  n 


as  — 9 **  and  \ - [C^  N ^] , where  we  have  set  a 1. 


Since  CL  t/  (Ct  nl) 7 


1 we  have 


(C.  Nx)  ft  (ai,  N.)  - g(ai)l— ^ n/o,  n(a.)),  i = 1,  1 


or 


Nr  [^(ai>  Ni)  ■'  N’(°.  n(a±)  / c^2). 


i = 1 , . . . , n . 


Since  rfa^,  tv  ) , i = 1,  ....  n are  independent,  we  have,  by  Theorem  3.2 
of  [1],  convergence  of  the  vector  process 


[r(a1,  N1)  - g/a^] 


N /2  ^a„»  Nn^  " g(jn)] 
1 n n n 


N(0,  n (31)/c^/2  ) 


M 


N(0,  n (a  )/C  z ) 
n n 


where  the  limit  random  v iriables  are  independent.  Using  the  Continuous 
Mapping  Theorem,  cf.  [1]  Theorem  5.1,  we  have 


Y(A) 


N1^2[?[a1,  Nx)  - g(ar)l 


VIf‘v  »„>  - 8<V 


N(0,  (£  **  <»>  n2(a[)/C1) 


1/2. 


i=l 


or 


N^[r(A,  Nr  ....  Nn)  - g(X)_]/[J  f j (Mn^a^/Cj72^  N(0,  1) , 


i“l 


Since  [C^  N^l/C^  — > 1 and  (a^,  N^)  —>  rj(a^)  in  probability, 


this  Imp  1 ics 


n 

[f>(A,  Nr  ....  Nn)  - g(A)J/i£  (A)  N,)/[Ci  fg]’/2  => 


N(0,  1) 


-5- 


[r  (A,  N1(  ....  Nn.)  - g(X)]/v-  (A,  j^,  ....  N^)  4 N(0,  U 
which  is  the  desired  result. 

The  above  proposition  allows  one  to  obtain  an  approximate?  100(1  -Y  )% 
confidence  interval  for  g(A)  as  follows.  For  sufficiently  large, 

i = 1 , ...»  n , 

P f-4>_1(l  - Y/2)  1 [r(A,  N , ....  N ),  - g(A)}/v(X,  , ....  N ) < 

(,  In  1 n — 


4>"1(1  - Y/2)  j » 1 - Y . 


This  may  be  rewritten  as 


P £ *(A,  Nr  ....  Nn)  - 0(X,  Nx,  ....  Nn)  $_1(1  - Y/2)  < g(A) 
r(A,  Nx Nn)  + v (A  , Nn>  - y./2)J  =l-y 


giving  the  desired  confidence  interval. 

Note  that  the  confilence  interval  obtained  at  A = 


= a^  reduces  tc 


r(at,  Nt)  + N~  2 fi(a1,  <J>_i(l  - Y/2)  , 


which  is  exactly  the  same  confidence  interval  which  could  be  obtained  based 
cn  the  M observations  at  A 5 a^.  Letting  Z (A  , N^,  N^)  denote  the 

length  of  the  confidence  interval  obtained  at  A , we  see  that 


UX,  N1?  .. 


, N ) - 2 v( A,  N. N ) 4*  (1  - y/2) 

n In 


( A),  l 2 (t 


-6- 


where  i (a,,  N.1  is  the  length  obtained  at  A = a.  . 

11  ' 1 


3 . Confidence  Band  Estimation  for  the  Function  g(A) 

As  noted  in  Section  2,  we  may  obtain  approximate  100(1  - Y )%  confidence 

. . l 

12i 

intervals  [A.,  A.]  for  g(A)  a".  A = a.,  i = 1,  ....  n: 


(1)  A^  <_  gCa^)  £ A2  with  probability  (1  - Y ) , 1=1,  n, 

where  A.1  = r (a, , N ) - n'  ^ n (a  , N ) $"1(1  - y /2)  and 

1 111  11  -L  , 


-l 


= fCa^,  N^)  + 2h(ai>  N^)  $ (1  - yjl).  Proposition  1 

enables  one  to  obtain  . onfidence  intervals  for  g(A)  at  any  point  ao-‘-Vr 

In  tit  section,  we  ohow  how  (1)  may  be  used  to  obtain  a confidence  band 

fir  the  entire  function  g(A)  over  a^  ^ A ^ an+i  • 

Before  stating  the  result,  we  .shall  need  some  additional  notation. 

1 2 

Define  the  column  vectors  A (m)  and  A (m) , m = 0,  1,  , , n as  follows. 

t h 1 2 

For  1 £ j n,  the  j components  of  A (m)  and  A (m)  are,  respectively. 


[• 

if 

j 

< 

Jl 

and 

m-j 

is 

evei. 

k 

i ,1 

j 

or 

j 

> 

m 

and 

j-m 

is 

odd. 

i ~i(m)  = 

i ■ 

< 

A? 

if 

j 

<_ 

ra 

and 

m-j 

is 

odd 

j 

1 

j 

or 

j 

> 

in 

and 

j-m 

is 

even 

i\  : ■ 

t 

l 

f 

A2 

if 

.1 

<_ 

m 

and 

m-j 

is 

even 

’ ■ o 

J 

or 

3 

> 

m 

and 

j-m 

is 

odd , 

~ j <tn)  " 

A* 

if 

j 

< 

m 

and 

m-j 

is 

odd 

j 

or 

j 

> 

m 

and 

j-m 

is 

even 

The  following  proposition  gives  the  desired  confidence  band  for  the  function 
g(A).  With  the  stated  probability,  the  inequalities  in  the  proposition  hold 
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xfan+1>  given  that  (1)  holds. 


simultaneously  for  all  values  a 


0 


PROPOSITION  2.  With  probability  (1  - Y^)  (1  - Y^  . . . (1  - Y ) , 


(2) 

for  a : A < a 


Y(A)  ^(m)  < g(A)  £ T ( A ) A2(m) 


tn  — 


m+1 


and 


m = 0,  1, 


. ,n. 


Proof.  Since  the  experiments  at  A = a , i = 1,  . . . , n are  independent, 
all  the  inequalities  (1)  hold  simultaneously  with  probability 
(1  - Yj)  (1  - Y2'  ...  (1  - Y ) . It  is  easy  to  see  that  the  inequalities 
(2)  imply  the  inequalities  (1),  since  f(a  ) Ak(m)  = Ak  for  k * 1,  2 and 
1 1 m 1 n.  It  is  therefore  sufficient  to  show  that  (1)  imply  (2). 

This  fact,  however,  is  proved  in  an  earlier  result  by  the  author;  see  [3]. 


The  above  proposition  may  be  interpreted  as  follows.  For  1 £ m £ n - 1 

th 

and  a < X < a the  bounds  for  g(A)  are  obtained  by  forming  two  (n-1) 
m nr*-i 

order  polynomials  passing  through  subsets  of  the  points  A^,  Ak  and 


4 


, A ' . The  upper  bound  polynomial  is  the  unique  polynomial  passing 
n 


2 2k 

through  A , A , , , and  A,,  1 < j < m-rl,  m + 2 < j < n,  where  k alternates 
m ra+l  j 

between  1 and  2 for  each  integer  step  away  from  the  interval  [m,  m +1], 


starting  with  k = 1 for  j = m - 1 or  j = ui  + 2.  Similarly,  the  lower 

bound  polynomial  is  the  unique  polynomial  passing  through  A^,  A^.,  and 

m nrri 

alternating  values  for  A ^ , l£j£m-l,  m+2£j£n,  starting  with  k ~ 

for  j ■ m - 1 or  j=m+2.  For  a^  A £ a^,  the  upper  bound  is  the 

2 12  1 

polynomial  passing  through  A^,  A^,  A^>  A^,  etc.,  whereas  the  lower  bound 

12  12 

polynomial  passes  through  A , A„,  A,,  A , etc.  Similarly,  for  a £ A £ 

1 l.  i 4 n 

21  2 1 

a^+^,  the  upper  bound  polynomial  passes  through  A^,  A^  An_2»  An_3’  etc*« 
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etc'. 


12  1 2 

whereas  the  lower  bound  polvnomial  passes  through  A , A A „ A ,, 

. n n-1  n-2,  n-J 


4 . Applications 

The  most  straightforward  application  of  Propositions  1 and  2 is  the 

situation  where  the  N observations  at  X = a^  are  independent  and 

identically  distributed  sample  variates  X.,,  ....  X..t  c , „ . 

J r il  iN^  from  a population 

with  mean  g(a^)  and  standard  deviation  n(a^).  Letting 

x N. 

r(a  N)=N JjX 
j-1 

and 

r i 1/2 

S(v  V ■ [sfri  ,£  % - *<v  \»2J 

we  know  from  the  Central  Limit  Theorem  and  the  Law  of  Large  Numbers  that 


Ni  - 8(ai)J/ri(ai)— ► N(0,  1) 


ri(a^,  Na)  — > nC^)  in  probability  , 

as  N — * cof  thus  satisfying  the  hypotheses  in  Proposition  1.  Similarly, 
Proposition  2 may  be  applied  once  confidence  intervals  are  obtained  at 
g(a^) , i ” 1*  • • • » n* 

We  shall  devote  the  remainder  of  this  section  to  a second  application 
which  is  based  on  the  regenerative  approach  for  analyzing  computer  simula- 
tion experiments.  (References  [4,  5,  6,  7,  10}  may  be  consulted  for  more 
details  on  the  regenerative  method.)  A basic  statistical  problem  in  simula- 
tions is  to  estimate  the  quantity  g(X)  = R ^f(X(X))^  where  f is  a general 
real-valued  function,  X is  an  input  parameter,  and  X(X)  is  the  stationary 
random  vector  associated  with  £x(s,X  ):s  > oj  , the  process  being  simulated. 
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0 * 


in 


In  the  regenerative  method,  we  observe  the  process  (X(s,  A);s  > ol 

random  cycles  of  lengths  (X),  a„(A),  — , a (X)  and  record  in  each  cycle 

' I Z N 

the  values  Y, (A)  Y„(A),  . ..,  i (A)  where  Y.  (A)  is  the  area  under  the 

1 Z N <c 

curve  f(X(s,A  ))  in  the  k cycle.  The  crucial  conditions  required  for  the 
regenerative  method  to  be  used  are  that  the  N pairs  {(Y  (A),a  (A)) 

K-  iC 

k = A.  2,  . . . , N}  are  independent  and  identically  distributed 
and  that  g(A)  = E {Y^(A)}  /E{a  ^ (A) } . These  conditions  hold,  for  example, 
in  simulating  a GI/G/S  queue  with  traffic  intensity  less  than  1,  letting 
cycles  be  defined  in  terms  of  the  beginning  of  successive  busy  periods. 

Now  define  the  column  vector  Uk(A)  = (Y^CA),  a^(A))  and  let 


1(A)  = 


X °tll(A)  ai2(A)  ' 


[a12 (A)  a22(A) 


denote  the  covariance  matrix  for  U^(A).  Denote  the  sample  mean  by 


U(A,  N) 


/ _ \ 

Y (A  , N)  A 


a(A,  N) 


s ZV) 

k=l 


and  the  sample  covariance  by 


i(A,  N) 


S11(A’  N)  S12(Af  N) 


s21(A,  N)  s22(A,  N) 


. N 

’ nTT  Z r\(A)  - U(A,  N)HUk(A)-U(A,  N)] 
k"l  " 


where  the  prime  denotes  transpose.  Next,  define  point  estimates  for  g(A) 
as  follows: 
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I — ~i 


( 1 ) Classical  est  1 mat o r 


f (X,  N)  - S , 
N) 


(2)  Beale  estimator 


rb(A,  N) 


= Y (a  , N) 
a(A,  N) 


[1  + s^(A,  N)/N  Y (A , N)  a(A  , N)] 

-L  -x-  

[1  + s22(A,  N)/N  a2 (A , N)] 


(3)  Tin  estimator 


f t (A,  N) 


Y(A,  N) 
a(A,  N) 


1 + 


s12(A,  N) 


s22(a,  N) 


Y(A,  N)  ct(A,  N)  a (A,  N) 


(4)  Jacknif e estimator 


fj(A’N)  = | kE  VA’ 

where  0 (A,  ,’)  = N[Y(A,  N)/a(A,  N>]  - (N  - 1)  [ £ Y (A)/  £ a (A)  ] 

j*k  J j?k  J 


Finally,  define 


n (A,  N)  = [s  (A  N)  - 2?  (A,  N)  s.. 
c 11  c 12 


(N)  + ?^(A,  N)  s22(A,  N)]/5(A,  N) 


and 


Hj(A,  N)  - Li  [0k(A,  N)  L fj(A,  N)]2/(N-1) j 


Now  let 


VA)  ■ VA) 


and  note  that  E{  Zk(A)}  - 0 and  define 
the  vectors  ^^(A),  k >_  1 j are  i.i.d.. 
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1/2 


- g(A)  A) 

o2(A)  - vat  | Zk(A)  | . Since 
it  follows  that  ^,Zk(A)f  k >.  1^ 


are  i.i.d.  By  the  central  limit  theorem  for  partial  sums  of  i.i.d.  random 
variables,  it  follows  that 


N 


£ Z (X)/Nl/2  o(X)  =4  N(0,  1) 


k=l 


as  N — which  pay  be  rewritten 

,1/2, 


N [f  (X,  N)  - g(X)1  ct(X,  N)/  a(X)  =~4n(0,  1) 
c 


Since  a(  X,  N)  — ^ E £ ot ^ (' ) } a.e.,  it  follows  that 


i /o  . 

N ' [rc(X,  N)  - g(X)]/n(X)  =4  N(0,  1)  , 


where  n(X)  = a(X)/E(a^  (X)}.  Now  it  may  be  shown  that 


N1/2[*  (X,  N).  - r,  (X,  N) ) — > 0 a.e., 

c b 


NV2[rc(X.,N)  - ft(X,  N)J 


0 a.e.. 


1/?  I ' , 

Mw*[r.  (X,  N)  - t (X , N)  ] — > 0 a.e., 

c J 

and  that  n (X,  N)  >n  ( X)  and  n,(X,N)  — > n(X)  in  probability  as 

C J 

N — > oo  . Hence*  the  conditions!  of  Proposition  1 are  satisfied  with  any 
of  the  following  substitutions: 


'cv  i’  ‘V 


■jx  i*  r 


To ■ illustrate  the  applicatio 


VI 

. w 

V 

or 

Vv 

V 

sub 

stituted 

for 

r(a 

. Mt). 

i * 1, 

2*  n i 

sub 

stituted 

for 

, n(ai, 

i - 1, 

2.  ....  n. 

h of  Propisitions  1 and  2,  consider  a 
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simulation  of 


the  customer  waiting  time  process  { W , n > 1 } in  a fil/C/l 

n — 

queue.  Suppose  we  wish  to  study  the  sensitivity  of  the  mean  stationary  wait- 
ing time  3(A)  = E { W ( \ ) } to  the  arrival  rate  A over  the  range 
2 1 A 1 6,  with  the  service  rate  U = 10.  In  what  follows,  we  shall 
illustrate  the  proposition  for  this  simulation  using  the  "classical"  esti- 
mators ?c  and  n though,  we  could  have  chosen  any  of  the  estimators  given 
abpve. 

In  the  queueing  simulation,  we  say  that  the  kth  busy  cycle  is  initiated 

with  the  arrival  of  the  customer  to  find  an  empty  queue.  Suppose  that 

simulation  runs  consisting  of  N = N = N = 10,000  busy  cycles  are  made  at 

I i.  3 

parameter  settings  A=ai=3,A=a  =4„  and  A = a = 5.  Let  a (A) 

L JR 

denote  the  number  of  customers  served  in  the  kth  busy  cycle,  and  let  Yk(A) 
be  the  sum  of  the  waiting  times  for  those  customers.  Then  it  may  be  shown, 
cf.  [4],  that  { (Yk(A),  ak(A),  k ^ 1 } are  independent  and  identically  dis- 
tributed, and  E(  W(A)}  * E{Y^(A)}  /E{ct ^ (a ) } • Hence,  it  is  appropriate 

to  apply  the  regenerative  method  as  discussed  above.  In  particular,  we  can 
compute,  for  i - 1,  2,  3,  Y(ai,N.),  aCa^  N^,  f^),  s^^,  N±) ,, 

and  s^^Ca^,  N^)  and  from  these  we  can  compute  rc(a^,  N^)  and.  nc  (c  ^ , N^) 

as  defined  avove. 

Suppose  that,  as  a result  of  these  computations  for  a simulation  run, 

r (a  N ) - .04  0 (a  » N.)  = .20 

c 1 1 c 1 1 

r c ( a2 » N2)  - .06  ( a2 • ly  - .30 

r (a,,  N ) - .10  n (a  N_)  - .50 

c 3 3 c 3 3 

Recall  that  approximate  100(1  - y)Z  confidence  intervals  at  A » a^  are 
given  by  y + 'Wa^  N^)  * _1(1  - Y/2)/N*/2  . For  example,  if. 
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95/o  confidence  intervals  are  desired,  these  wou'd  be 


.03608  _<  g ( 3 ) £ .04392  with  95%  confidence' 

.05412  1 g(4)  1 .06588  with  95%  confidence 

.09020  £ g(5)  £ .10980  with  95%  confidence. 

Now  assume  that  g(A)  = E (W(A)}  is  approximately  quadratic  for  2 ± A <_  6. 
From  Proposition  1,  an  approximate  100(1  ~Y)%  confidence  interval  at  A is 
given  by  9(X,  N^  N3>  + G(X,  l^,  Nj)  *_1(1  - Y/2),  2 £ x £ 6, 

where 


r(A,  Nx,  N2,  N3) 


r.  04 


*(■*> 


.06 


.10 


v(X,  Nr  N2,  N3)  = J(X)(.20)2  + Y 2 (A)  (.30) 2 +V  2(A)(.50)2J  1/2/100 

and 


YfA)  * (1  X A2) 

13  9 

1 4 16 

-l 

= (1  X A2) 

10  -15  6 

-9/2  8 -7/2 

1 5 25 

L 

1/2  -1  1/2 

Thus,  for  example,  a 95%  confidence  interval  for  g(3.5)  is 

.04269  ± g(3.5)  ± .05231  with  95%  confidence. 

Finally,  using.  Proposition  2 together  with  the  above  confidence  intervals  at 

3 

A » 3,  4,  and  5,  we  obtain  a 100(.95)  % confidence  band  for  the  function  g(A) 
over  2 £ X £ 6; 


for  a < \ 

m ■ 


Y(A)  A1(m)  <^g(A)  <_  Y (\)A2(m) 

< a and  m - 0,  1,  2,  3,  wh*»re 
— nrt-1 
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A^O)  = 


Ai(2)  * 


.03608 

.04392 

.03608 

.04392  | 

.06588 

a2(0)  = 

.05412 

A 1 ( 1 ) = 

.05412 

a2(D  = 

.06588  | 

.09020 

.10980 

- 

. 10980 

. 09020 J 

— — 1 

.04392 

.03608 

r ~i 

.03608 

j 

~ ! 

.04392 

.05412 

.09020 

A* (2)  = 

.06588 

.10980 

Al(3)  = 

.06588 
. 09020 

a2(3)  = 

1 

.05412 

.10980 
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